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STEADY-STATE SOLUTIONS IN AN ALGEBRA OF GENERALIZED 

FUNCTIONS: 

LIGHTNING, LIGHTNING RODS AND SUPERGONDUCTIVITY 

TODOR D. TODOROV 


Dedicated to the memory of Christo Ya. Christov* * on the 100-th anniversary of his birth and 
to James Vickers on the occasion his 60-th birthday. 

Abstract 

Formulas for the solutions of initial value problems for ordinary differential equations with 
singular J^^^-like driving terms are derived in the framework of an algebra of generalized functions 
(of Colombeau type) over a field of generalized scalars. Some of the solutions might have physical 
meaning - such as of the electrical current after lightning or under superconductivity - but do not 
have counterparts in the theory of Schwartz distributions. What is somewhat unusual (compared 
with other similar works) is the involvement of infinitely large constants, such as <5(0), in some 
of the formulas for the solutions. 

1. Introduction and Notation 

1. Our notations are similar to those in Bremermann [1] and Vladimirov m- N,No,M, M+ 
and C denote the set of the natural, whole, real, positive real and complex numbers, 
respectively. f(M) = C°°(M) stands for the class of C°°-functions from M to C and 
D(M) denotes the set of test-functions on M. We denote by D'(M) the space of Schwartz 
distributions on E and by P(|_(E) the space of all distributions in D'(E) supported by 
the interval [0, oo). 

2. The main purpose of this article is to derive formulas for the solutions of the initial value 
problem: 

(1) P{d/dt)y = f{t), y{0) =y'{0) = ■■■ =y^^~^\0) =0, 

in an algebra T(E) of generalized functions over a field of generalized scalars C (of 
Colombeau type) constructed in Todorov & Vernaeve [13]. Here P G C[x] (and more 
generally, P G C[x]) is a polynomial in one variable of degree k and P{d/dt) stands for 
the corresponding differential operator with constant coefficients. The right hand side of 
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the equation, /, is a distribution (and more generally, a generalized function in £1(M)). 
The solution of the intitial value problem ([1]) are often called steady-state solution of 
the differential equation P{d/dt)y = /(t), which explains the title of the article. We 
are mostly interested in initial value problems ([T]) which do not admit distributional 
solutions even when P E C[x] and / is a distribution with compact support. 

3. Recall that if T E P'(M) is a distribution and xq is a real number, then the value 

T(xo) E C exists in the sense of Robinson ([ID], §5.3) if and only if T is a continuous 
function on a neighborhood of xq (Vernaeve & Vindas [E]). Consequently, ([I|) admits a 
(unique) distributional solution y if and only if the equation P{d/dt)y = f{t) admits a 
(particular) solution yp{t) E for some e E M_|_. Here is a typical example of 

differential equation which does not admit a steady-state distributional solution. 

Example 1. The IVP y" + uj‘^y = 8'{t), y(0) = y'(0) = 0, where oj E M_|_, does not admit 
a distributional solution. Indeed, the general solution of the equation y" -\- uf^y = 6'{t) 
is presented by the family y = ci cos cut -|- C 2 sina;t -|- Gg(t), ci,C 2 E C, where Go{t) = 
l^H(t) sin ujt is the Green function of the operator ^ . Here H E P'(M) stands 

for the Heaviside distribution with a kernel, h{t), the unit step function. We also have 
GQ(t) = H{t)cosu)t and ^(((t) = 5{t) — w^Go(t). The initial condition y{Q) = 0 implies 
Cl = —Gq(0) and y'{0) = 0 implies C 2 = —^[5(0) — a;^Go(0)] = —since Go(0) = 0. 
Thus the only candidate for a solution of our IVP is y = —Gq(0) cos cut — ;^<5(0) sinujt 
H{t) cos ujt. However, neither 6(0), nor Gq{ 0) exists in the theory of distributions. 

4. Strangely enough, it seems the above phenomenon - non-existence of a distributional 
solution - is widely unknown (if known at all); way too often initial value problems for 
steady state solutions without distributional solutions (as those presented above) appear 
in the mathematical literature. We have selected three typical examples (among many): 
(a) In (Edwards &: Penney |6|, 7.6 Problems) the problems # 1, 2, 4, 8 do not admit 
distributional solutions and are supplied with wrong answers. In (Schiff |9|, p. 104) the 
function g{t) = sinht is presented as the solution of g" — g = 5{t),g{t)) = ^'(0) = 0 
(although ^'(0) = 1) and in the recent article (Strang [T2|, p. 1245) the function g{t) = 

is presented as the solution of g' — ag = 5{t), g{0) = 0. The authors obviously 
are unaware that these initial value problems do not admit distributional solutions (and 
never bothered to verify the correctness of the final result). 

5. It is very tempting to try to solve the initial value problem ([T|) by the Laplace trans¬ 
form method and this is what the authors cited above have tried to do. We should 
mention that: (a) The Laplace transform will always produce the correct solution to ([I|) 
if dlj) has a distributional solution, (b) If ([1|) does not admit a distributional solution, 
the Laplace transform method (if applied correctly) will produce a particular solution 
of the differential equation which does not satisfy all initial conditions. However, the 
Laplace transform will not alarm us for non-existence of solutions (verifying the final 
result is always advisable), (c) Strictly speaking the formulas in the usual tables of 
Laplace transforms (which appear everywhere in the textbooks and also in Wikipedia) 
are logically inconsistent (Todorov [H]). For example, these tables suggest wrongly that 

= sint instead of the correct one = H{t)sint. For the readers 

who prefer to use the Laplace transform (over the Fourier transform) for the purpose of 
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finding a particular solution of a differential equation, we recommend the table ’’Laplac- 
eTable” (which incudes as well) which can be found in the webpage of the author 
of this article. Alternatively, the reader might refer to to the inversion integral formula 
for (Folland [1]), §8.2). 

6. Notice that the above formula y = —G'fy{0) cosojt — ^|;^sina;f + H{t) cos ujt (although 
meaningless within the Schwartz theory of distributions) make perfect sense in any of the 
numerous algebras of generalized functions of Colombeau type (Colombeau [3]). This is 
because the value T(0) is well defined - as a generalized number - for any distribution 
T. The main goal of this article is to show that these formulas (and many similar to 
them) not only make sense, but they are actually the solution (steady-state solution) of 
the initial value problem from which these formulas originate. 

7. Except the example at the end of our article, the solutions of ([T]) in this article are linear 
combinations of Schwartz distributions with coefficients in M. The fact that the set of 
scalars C of the algebra of generalized functions T(M) (Todorov &: Vernaeve [13]) is a 
field (an algebraically closed non-Archimedean field) - not a ring with zero divisors as in 
Colombeau m) - is important for our approach; it allows us to transfer the familiar linear 
algebra arguments from the classical theory of ordinary linear differential equations with 
constant coefficients to the similar equations in £1(M). 

8. In the last section of the article we endow some of our initial value problems with phys¬ 
ical interpretation inherited from the Kirchoff law: LI" R I' + ^ I = V'{t), 1(0) = 

0, LI'{0) = 0, for an electrical Li^C-circuit with inductance L, resistance R, capac¬ 

itance C, driving electromotive force V{t) and (steady-state) current I{t) (Bremer- 
mann [I], Ch.lO). We derive formulas for I{t) for some extreme and violent physical 
phenomena such as lightning and superconductivity which do not have counterparts in 
classical analysis and the Schwartz theory of distributions. We challenge the common 
wisdom that only the real-valued functions have the right to present physical quantities: 
our formulas for the electrical current I{t) are often generalized functions with infinitely 
large values which sometimes keep their infinitely large values in time. 


2. Linear Independence in T(M) 


We ask the reader to refer to (Todorov & Vernaeve [13]) for the construction of the algebra 
^(R) and the field of its scalars C. Here is a summary of this construction: 

1. If S' is a set, then every function of the form / : P(R) —>■ S will be called a (/?-net in S 
and will be denoted often by (ftp), where cp G P(R) is treated as a parameter (Todorov 

&: Vernaeve [13], §3). The elements of S(R) are equivalence classes q{fp) (denoted also 
by fip) of (/9-nets in S(R). Similarly, the elements of C are equivalence classes q{Ap) 
(denoted also by A^) of (^-nets in C (Todorov & Vernaeve [l3], §4). In particular, the 
generalized number s = q{Rp) is called canonical infinitesimal, where 


( 2 ) 


sup{||x|| : X G R, (p{x) 0}, 

1 , 


V 7^ 0, 

(^ = 0 . 


stands for the the radius of support of (/? G P(R). 
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2. The set of the (generalized) scalars C = {/ E £1(M) : /' = 0} of £1(M) is a field and we 
also have C = M(z), where M is the field of the real (generalized) scalars of £1(M). Both 
C and M are non-Archimedean fields; C is an algebraically closed field and M is a real 
closed field (Todorov Vernaeve [l3]). If x E C, then x is infinitesimal if |x| < 1/n for 
all n E N. In particular, s is a positive infinitesimal, i.e. 0 < s < 1/n for all n E N. 
If X, ?/ E C, then we write x ~ 7/ if x — ?/ is infinitesimal. We denote by the set (an 
integral domain) of the finite elements x E K for which |x| < n for some n E N. 

3. We denote by l : P'(M) —)• £1(M), where i{T) = q{T the canonical embedding (of 
Colombeau type, Colombeau m) of the space of Schwartz distributions into f(M) (we 
shall use here i instead of used in Todorov &: Vernaeve m, p. 222). Here T * (/j 
stands for the convolution product between T and cp in the sense of theory of distributions 
(Vladimirov [T7j, 79-80). In particular, i{6) = q{p>) (Todorov & Vernaeve [13], §5). 

4. The elements / of £1(M) are pointwise functions of the type / : —>■ C. The latter holds, 

in particular, for / = i{T), where T E P'(M). If / E £^(M), then i(/) is an extension 
of / in the sense that i(/)(x) = /(x) for all x E K. On this ground we shall often use 
the notation / instead of the more precise i(/) leaving the reader to figure out what 
we actually mean from context. For each A; E Nq the value i((5^^^^)(0) is infinitely large 
number in M, i.e. n < |i((5^^^^)(0)| for all n E N. Also, i((5^^^'''^^)(0) = 0 (since t{6){t) is 
an even function). 

5. The algebra £1(M) and its scalars C admit also an axiomatic definition (Todorov [T3| and 

Todorov |15] ). We should mention that under some assumption £ (R) is isomorphic to the 
algebra of p-asymptotic functions p£(JR), defined in (Oberguggenberger &: Todorov |8|) 
and R is isomorphic to Robinson’s field of p-asymptotic numbers ^R (Robinson [llj . 
Lightstone & Robinson m)- We prefer to work with £1(R) (rather than with ^£1(R)), 

only because the embedding l : P'(R) —)■ f(R) is canonical (while the similar embedding 
in ^£1(R) depends on a fixed non-standard mollifier). 

Definition 1 (Vector Spaces span(S')). (1) For every S C P'(R), we let spah(S') = 

{ I'iTn) : m E N, E 5, Ori E C}, the span of l{S) within the algebra £ (R). We 

supply spah(S') with the operations of a differential vector space over the field C inherited 
from the algebra ^(R). 

(2) We supply spah(P()_(R)) with the convolution product * : P(^_(R) x si)ah(P(|_(R)) —>■ 
spah(P()_(R)) by T -k f = i^iT kTn), where T -kTn G VfifR) is the convolution 

product between T and Tn in the sense of the theory of distributions (Vladimirov [T7], p. 

77). 

We observe that * f = for n = 0,1,2,... . Notice that for the value (T * /)(0) E C 
we have (T* /)(0) = ^(F*r„)(0). Recall that l{T *T„)(0) = q{{TkTnk(p)(0 )), where 

q[{T kTnk <p)(0)) is the equivalence class of the <p-net ((T kTnk (p)(0)) in C (Todorov & Ver¬ 
naeve m, 213). 

In what follows W{fi,..., fk) denotes the Wronskian of the functions /i,..., /^ and we denote 
by IF(/i,..., /fe)(x) the value of IF(/i,..., fk) at x. 
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Lemma 2 (Linear Independence in ^ (R)). Let yi,... ,yk E £^(R) be solutions of the homogenous 
equation P{d/dt)y = 0 for some polynomial P G C[x] (in one variable with complex coefficients) 
of degree deg(P) = k. Then yi,... ,yk are linearly independent in the C-vector space f(R) if 
and only if their i-images, i{yi ),..., i{yk), are linearly independent in the C-vector space f’(R). 

Proof. i{yi ),..., i{yk) are solutions of P{d/dt)y = 0 in ^(R) since i is a differential ring embed¬ 
ding of £1(R) into £1(R) and i{W{yi,.. .,yk)) = W{L{yi),.. .,L{yk)) (for the same reasons). Also, 
W{yi,..., yk){x) = W{i{yi ),..., i'{yk))ix) for all x G R since each i{yn^) is a pointwise extension 
of yn\ Now, yi,... ,yk are linear independent in £1(R) if and only if W{yi,... ,yk){0) 7 ^ 0 if 
and only if W{L{yi),... ,i{yk))iO) 7 ^ 0. The latter is equivalent to the linear independence of 
i(yi),..., L{yk) in T(R) by the usual arguments since (ci(y„))W( 0 ) = cyn\o) for all c G C and 
all n and i. □ 


3. Associated Generalized Constants 

Definition 3 (Associated Generalized Constants). Let P G C[x] be a polynomial in one variable 
of degree k and let Gp G P'(R) he the Green function of the operator P{d/dt) : P'(R) —)■ P'(R), 
i.e. the distributional solution of the initial value problem P{d/dt)Gp{t) = (I(t), Gp{t) = 
0 on (— 00 , 0 ). Let i{G^p'^) be the i-image of G^p^ into T(R). We say that i(Gp^)(0) G C, 
n = 0,1, 2,... , /c — 1, are the generalized constants associated with the operator P{d/dt). 


Here are three examples - presented as lemmas - of particular differential operators along with 
their generalized constants. 

Lemma 4 (Green Function of + b). Let a ,6 G R, o 7 ^ 0 and let g G P'(R) be the Green 
function of the operator + b, i.e. ag' + bg = d{t), g{t) = 0 on (—oo,0). Then: 

(i) : g{t) = ^H{t)e~a^ , , where H G P'(R) is the Heaviside distribution with kernel h, the 
unit step-function. Consequently, g'{t) = ^(d{t) — bg{t)^. More generally, g^'^\t) = 

\ ait), n = i, 2,^... 

(ii) : The only associated generalized constant is t(fl')(0) G R and t.{g){0) ~ l/2a. 


Proof, (i) The general solution of ay' + by = 6(1) in P'(R) is given by the family y{t,c) = 
gg-W/a _j_ ^ where c G C. The initial condition y(t) = 0 on (—oo,0) implies c = 0. 

Thus g{t) = ^Fr(t)e“a*. The rest of the formulas follow by direct differentiation. 

(ii) We have i(ff)(0) = qig*‘p)i0)) = ^q{ (^(x) dx) = lq{ (p(x) dx) (Todorov& 

Vernaeve m)- On the other hand, by the mean value theorem for integration (Hobson [3]) for ev¬ 


ery test function ip there exists in {—Rp, 0 ) such that e ^ <p(x) dx = e ffp t(x) dx+ 

!x^Tix)dx. Thus/,( 5 )( 0 ) = \[qie-^)qijplpix)dx)+qiJl^ipix)dx)] « 1 g((^(x) dx + 
+ !x^Tix)dx) = ^q(f^^(p(x)dx) = ^ since q{R^p) = s is a positive infinitesimal, thus 




qie 


bliRip) 


= e 


= e 


1, amd qi^jf'p ifix) dx) is a finite number. 


□ 
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Lemma 5 (Green Function of + b-^ + c). Let a, 6, c G M, a / 0 and 6^ — 4ac < 0 and 
let G G be the Green function of the operator + b-^ + c, i.e. aG" + bG' + cG = 

6{t), G{t) = 0 on (—00,0). Then: 

(i): G{t) = ^ sin (cut), where a = b/2a and uj = ^y/Aac — b'^. Consequently, 

G'(t) = -^H(t) \ — Q;e““* sinuit + oj e““* cos cut] = —aGit) + — H(t)e~°‘^ cos cut, 
aui a 

G"{t) = - 5{t) — - G'{t) — - G{t) = - 6{t) — X: H{t) e““* cos cut + — -^ H{t) e“"* sin cut]. 

a a a a aoj 


(ii): Both generalized constants t(G)(0) and l{G'){ 0) are generalized numbers inM such that 
t(G)(0) 0 and i(G')(0) l/2a. More precisely, |i(G)(0)| < s/2a and |t(G')(0) — ^| < 

where s = q{R^p) is the canonical infinitesimal o/M. 


Proof (i) We refer the reader to Bremermann [T], p. 127. 

(ii) We have l{G) = q{G-k ip) = smoj{t — x) (p{x) dx); thus we have 

i(G)(0) = q{ — flji sin (uix) (p{x) dx). On the other hand, by the mean value theo- 


, there exists x^p G {—R,p, 0) such that — ^ e“* sin (cox) p(x) dx = 


rem 


_ _\_p—aR, 

auj^ 

q-ols sin(a;s) 


sin(wi?,p) p(x)dx. Thus |t(G)(0)| < (/(I 


^-aRip 


|/!L < 




^ sin(a;s) I rO 




’ sin(aj s) 1 _ e 

aco 2 


^ sin(a;s) 


2a u 


< To 


show that t(G')(0) l/2a, it suffices to show that i(H(t)e “*cosa;t)(0) « 1/2. Indeed, 

i(H(t)e~°‘^ cosut) = q((H(t)e~'^^ cosujt)-k p) = q(^ cosu;{t — x)p(x) dx). Thus 

i( F7(t)e““* coscut)(0) = g( coswx (/?(x) dx) =q(f^j^ 6°"^ coscox p(x) dx). On the other 

hand, by the mean value theorem for integration (Hobson [5]), there exists x^ G (—i7<^,0) such 
that e°^^ cos LUX p(x) dx = e~°^^r cos (cuRp) p(x) dx + p(x)dx. Thus 
i(H(t)e~°‘^ cos out) (0) = q(e~°‘^r cos cuRtp)q( p{x) dx)+q(^ p(x) dx) ^ g( p(x)dx)+ 
q( p(x) dx) = q{ t(x} dx) =1/2 since q(e~°‘^r cos (cuRp)) = e““^ cos (cos) ps 1 . □ 


Lemma 6 (Green Function of ^ Tcu^). Let cu G M+ and let Gq G P((_(M) be the Green function 
of the operator ^ +a;^, i.e. G'/ + oj‘^Gq = d(t), Go{t) = 0 on (—oo,0). Then: 

(i): Go{t) = ^H(t) sin (cut), where (as before) H G P'(M) is the Heaviside distribution 
with kernel h, the unit step-function. Consequently, G'^ft) = H(t) cos (cut), GQ{t) = 
6{t) —uP'Go^t) and ^(('(t) = 5'(t) — uj‘^GQ{t). More generally. 


^(^2n+2)(^) ^ ^^_^^k^2k^{2n-2k)^p^ ^ ^^„+i^2(n+l) 

k=0 

n 

^(^2n+3)(^) ^ ^ ^^„+i^2(n+l) 

k=0 


for n = 0,1,2,... . 
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(ii): The generalized constants i(Go)(0) and c(Gq)(0) are generalized numbers in M such 
that t(Go)(0) Ri 0 and i(GQ)(0) = 1/2. More precisely, l{Go){ 0) < sjl, where s = q{Rip) 
is the canonical infinitesimal o/M. Consequently, we have 

n 

^^^(^ 2 n+ 2 ))(o) ^ ^ ^2{n+l) ,(Gp)( 0 ), 

Proof. Except for the strict equality t(GQ)(0) = 1/2, the results follow from the previous lemma 
for a = 1,6 = 0 and c = We have i(GQ)(t) = q(f^^ cos uj(t — x)(p(x) dx); thus we have 
''(G'o)( 0 ) = q(f^^cos(cvx)(f(x)dx} = qQ cos (cvx) (p(x) dx) = i fid)(x) cos(u)x) dx = 

In the last two formulas we have taken into account that = 0 for odd n. □ 

4. Steady-State Solutions in £^(M) 

Theorem 7 (First Order ODE). Let a,b G R,a 0 and let g{t) = ^H(t)e~a^ be the Green 
function of the operator a^ + 6 (Lemma^. Then for every f G spaii(P(,_(]R)) (Definition\fi^ the 

initial value problem ay' + by = f{t), y{0) = 0 , has a unique (steady-state) solution y in £^(]R) 
m, given by the formula 

(3) hit) = {g* f){t) - {g-k f){0)e~^\ 

where gk f G spah(D((_(M)) and {g* f){0) G M. 

Proof. We have / = Yl^=i I'iTn) for some m G N, G V'^iR) and an G Chy assumption. 

Next, the general solution of the equation ay' + by = Tn in P'(M) is given by the family yn{t, c) = 

{gkTn){f) + where c G C. Consequently, the general solution of the equation ay' -\-by = 

L{Tn) in £^(M) is given by the family yn{t,c) = fig k Tn){t) -|- ce~a', where c G C. The initial 
condition y„(0, c) = 0 (treated as a equation for c in C) produces Cn = —fig k Tn){0). Thus 
the solution of the initial value problem ay' + by = fiTn), y( 0 ) = 0 , in T(M) is yn{t,Cn) = 
fig k Tn){t) — fig k r„)(0) e“ a*. The latter implies that y = Unit, Cn) is the solution of 

ay' + by = f{t), y{0) = 0 by linearity. Thus y = Y.n=i (^n{fig * Tn){t) - fig k T„)(0) e"«*) = 

(g* Dfi) ~ ( 5 */)( 0 )e~“* as required. □ 

Theorem 8 (Higher Order ODE). Let P G C[x] and deg(P) = k. Let f G spah(P()_(M)) 

(Definition Then the initial value problem ([I]) has a unique (steady-state) solution y in 

the algebra T(M) given by the formula y = X]n=i fiVn) + G k f, where yi,... ,yk are linearly 

independent solutions of the homogenous equation P{d/dt)y = 0 in T(K), G G D(,_(M) is the 
Green function of the operator P{d/dt), the convolution product Gk f G si)ah(D(|_(M)) is in the 
sense of Definition{^ and the generalized constants Cj G C are determined by 

W{yi,..., yn-uG k /, yn+i, ...,yk) (0) 

W{yi,... ,yn-i,yn,yn+i, • • • ,yfc)(0) 


(4) 


n = 1,..., A;. 
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Proof. We observe that for each n the family Yn{t, Ci,..., Ck) = CiUi + G -k Tn, where 

Ci G C, presents the general solution of the equation P{d/dt)y = Tn in P'(M). Consequently, 
the family yn{t, ci,..., c^) = <'{yi)P‘'{GkTn), where Ci G C, presents the general solution 

of the equation P{d/dt)y = r(T„) in £1(R) since i{yi ),..., i{yk) are solutions of the homogenous 
equation P{d/dt)y = 0 in £1(M) which are linearly independent in £1(M) by Lemma [2l Thus the 
family 7/(t, ci,..., Cfc) = Cn i{yn) + G * /, where c* G C, presents the general solution of 

the equation P{d/dt)y = G * / in T(R) by the linearity of the operator P{d/dt). The initial 
conditions y(0) = • • • = = 0 leads to the system of linear equations for the constants 

Cl,... ,Cfc in C: 

ELiCn2/n(0) =-(G*/)(0 ), 

< . 

since i{cn yn^)(0) = Cn yn\G) for each i and n. The Crammer rule produces the formulas dH for 
constants Cj in the the formula y = G k f + En=i 

Remark 1 (Simplified Notation). We shall often drop the l’s by letting yi = i{yi) on the ground 
that i[yi) are pointwise extensions of yt. In this simplified notation the above formula will be 
written simply as y = Yln=i CnUn + G-k f. 

5. LRC-Electrical Circuit: Lightning, Lightning Rods and Superconductivity 

We present several examples of steady-state solutions to first and second order ordinary 
differential equations with constant coefficients. We supply these initial value problems with 
physical interpretation borrowed from Kirchoff law for an electrical LRC-circle. Some reader 
might prefer to ignore physics and focus on mathematics. 

Definition 9 (Lightning, Lightning Rods and Superconductivity). Let V{t) G spah(P()_(M)) 
(Definition Let L,R,G G M such that L > 0,R > 0, G > 0 and Lf + R? ^ 0. Let 
I{t) G TW the solution of the IVP: LI" pRI' + ^I = V'{t), 1(0) = 0, L/'(0) = 0. Let 
A G R+. Then: 

1. We associate the above initial value problem with the Kirchoff law for an electrical 
LRG-circle with inductance L, resistance R, capacitance G, driving electromotive force 
(voltage) V{t) and (steady state) current I(t). 

2. V{t) = Ai{H){t) is called switch with amplitude A: 

3. Let n G No. Then V{t) = A L{dG)){t) is called lightning of order n with amplitude 

A. 

4. The operator L-^ R-^ + ^ is called a lightning rod if L,R and G are infinitesimals 

(G must be a non-zero infinitesimal). 

5. If R is infinitesimal (R = t) is not excluded), we talk about superconductivity. 
Lemma 10 (Steady-State Current). Let L,R,G and V be as in Definition\^ 

(i): If L = 0, then the Kirchoff law reduces to RP ^ I = V'{t), 1(0) = 0 and its 
solution in T(M) is given by the formula I{t) = {gkV'){t) — {g kV'){G) e~Rc ^ where 
9{t) = ^H{t)e~T^. 
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(ii): Let E? — ^ < Q. Then the Kirchoff law reduces to LI" + RI' + ^ I = V'{t)^ /(O) = 
I'{tS) = 0, and its solution in £^(M) is 


/(t) = —(G * y')(0) e 2 l* coswt— 


RiG*V')iO) + 2L{G'i^V'){0) _j? 


2Luj 


e 21 , sinwi-l- 


+ {G*V’)it), 


where oj = %■ — t'{G){0) ~ 0 and t{G'){0) ^ 1/2L (LemmalS^ and G{t) = 

-^H{t)e~^^ sinwi (consequently, G'{t) = —■^G{t) + ^ H{t) e~^^ cosLut, and G" = 

^5{t)-^G'{t)-^G{t) = ^5{t)-^ H{t)e-^^ cosujt+-uj'^) H{t) e-^^ sincat). 
(iii): Let R = 0 (superconductivity) and L,G & M+. Then the Kirchoff law reduces to 

L I" + ^ I = V'{t), 7(0) = I'{0) = 0, and its solution in £^(M) is given by: 

, (Go*'l7')(0) (G'oAy')(0) . 

I{t) = - - - - cos ivt --- - sin uit + — (Go * V')(t), 

1j IjUJ 1j 

where io = i'(Go)(0) ~ 0, i-(Gq)(0) = 1/2 (Lemma\^ and Go(i) = ^H{t) sincjt 

(consequently, G'^ff) = H(f) cos Lot, and Gg = 6{t) — Go(t) = 6{t) —LoH{t) sinoii/. 


Proof, (i) is identical to Theorem [7] for f(t) = V'(t). 

(ii) Vh(e“^*cosa;t, e“^*sina;t)(0) = lj, 1T((G ★ 17'), e“^* sina;t)(0) = a; (G a 17')(0) and 

W{e-^^cosLot, (G A 17'))(0) = ^(G* F')(0) + (G' a y')(0). Thus ci = -(G a V’)(0) and 
C 2 = -hi^iG* V'){0) + (G' A y')(0)) = _ ^(g*^')(0H2L(GGi/')(0) follows by 

Theorem [HI 

(iii) is a particular case of (ii) for R = 0. □ 


Example 2 (L = 0 and Switch). Let V{t) = Ai(H)(t) (switch), where A G M+. The corre¬ 
sponding initial value problem is RI' + ^ I = A i(6){t), 1(0) = 0. ITe have G-k f = g-k{A i{5)) = 
AL{g) and Lemma [TU produces the steady-state solution I(t) = A[L(g)(t) — t( 5 )( 0 )e“flc] or, 
equivalently, I{t) = A[^ L(H(t)e~Rc) — /,(g')(0) e“ , where L{g){0) is a generalized constant 
in M which is infinitely close to 1/27?. 


Example 3 (L = 0 and Lightning of Zero Order). Let V{t) = AL(6)(t) (lightning of zero 
order). The corresponding initial value problem is RI' + ^7 = Ai(5'){f), 7(0) = 0. ITe have 
Gk f = g k (Al(6')) = AL(g') = ^ [i(6)(t) — ^t{g){t)] and Lemma UTK vroduces the steady-state 
solution 


m = 71 


4<^)(0) 4g)(0) ,-4. 


where i(g)(Q) and t((l)(0) are generalized constants inU. such that L(g)(0) « 1/27? and i-(6)(0) is 
an infinitely large constant. 


Example 4 (L = 0 and Lightning of Zero Order, Superconductivity). Let a,b & M+. ILe 
are looking for the steady state current I(t) in an electrical LRG-circle with inductance L = 
0, infinitesimal resistance 7? = a/i((5)(0) (superconductivity) , infinitely large capacitance 
C = i{6){0)/b and driving electromotive force V(t) = i(6){t) (lightning of zero order), i.e. 













10 


TODOR D. TODOROV 


I' + ^ H^) — 0; or equivalently, a I' + bl = l{S){0) i{6'){t), 1(0) = 0. We 

apply the Lemma [lU for V = i((5)(0) i.(<5') and the result is g-kV' = (^t(6)(t) — be(g)(t)). 

Thus the steady state solution is 


Iit) = 


4<5)(0) 


— i(5)(0)e “* + 6i((j()(0) e - i(^H{t)e “*)+i(5)(t) 


Example 5 (L = 0 and Lightning of n-th Order). Let A G M+ and V{t) = At{6^^^){t) (light¬ 
ning of n-th order) . The corresponding initial value problem is RI'I = AL{6^'^~^^'^){t), 1(0) 
0. We have {Gk f){t) = A t{g^'^~^^^){t) and with the help of Lemma^we obtain: 


A 


(-1)^+1 r((5("-^))(0) (-l)’^r(ff)(0) 


k=0 


{RCf 


e Rc + 


+ ,^l^,A {H{t)e-Rc) + 


{RCY+^ 


k=0 


(-1)^ 

{Rcy 


RnCn+l 

i(j(”-^))(t) 


e RC + 


where (as before) L{g){0) « 1/2R and r((5^’^^)(0) are infinitely large numbers and i(5^’*'''^))(0) = 0 
/orn = 1,2,... . 

Example 6. Let ^ — R^ > 0 and let A G M+. 

1. Let V{t) = Ai{H){t) (switch). The solution of LI" + RI' + = Ai{5){t), 1(0) = 

I'iff) = 0 is given by 

I{t) = A —i{G){0)e 2i*cosa;t— (--—-)e 2 L*sina;t+ 

L \ oj 21 j(jJ / 

1 R ~\ 

+ -—r(Lr(t) e“2r* sinwt) . 

JjCJ 

2. Let V{t) = AL{5{t)) (lightning of zero order). The solution of LI" + RI' ^ I = 
Ai{6'){t), 1(0) = I'{0) = 0, is given by 

I{t) = A — i(G')(0) e“^*coswt+ 


, t(5)(0) Rl{G'){0) r(G)(0)\ r, . 


1 R R 

+ y r [H(t) e“ 2 i* cos uit — I'iRit) 

1j 21j^UJ 


e zi^sinw 


t)) 


3. Let V{t) = AL{d'){t) (lightning of first order). The solution of LI" + RI' + y I = 
Ai{6"){t), 7(0) = /'(O) = 0, is given by 


i^t) = a[[- ^ 


e 2L*cosu;f+ 

R , 

sin ojt+ 

l{H( t) e~ 2 r * sin ivt') 


L ' L ' LG J 
Ri{6m , {2L-R^G)i{G'm iir(G)(0)^ _r, 

2L^uj + 21^ ^ 

i{6){t) 




2L^Guj J 
R^C - 2L 
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where lo = ^ , ^(G)(0) « 0 and l{G'){0) ^ 1/2L (Lemma\^. Notice that 

l{S'){0) = 0 (Todorov & Vernaeve [T3] ). 

Example 7 (Superconductivity). Let i? = 0 (superconductivity), uj = Xj^/LC , A G M+ and 
Go{t) = ^H{t) sin (iot). Recall that /-(Go)(0) « 0 and i(GQ)(0) = 1/2 (Lemma\^. 

(1) IfV{t) = Ai{H){t) (switch), then the solution of L I" + ^ I = Ai{5){t), 1(0) = I^O) = 

0, is I{t) = ^ ^ sinwt) — ^ sinwt — t(Go)(0) coswt . 

(2) IfV{t) = Ai{6){t) (lightning of zero order), then the solution of L I" + ^ I = 

Ai{6'){t), 7(0) = I'{0) = 0, in 7’(M) is I{t) = ^ — sin ujt + t{H{t) cos cot) — 

^ cos cut + w i(Go)(0) sina;t . 

(3) If V{t) = Ai{5'){t) (lightning of first order), then the solution of L I" + ^ I = 

AL{6"){t), 1(0) = I'{0) = 0 in£{M) is I(t) = ^ cos u}t+^ sinuit—to l{H( t) sinujt)+ 

uj'^ l{Go){ 0) cos iot + L{6){t) . 

(4) IfV{t) = A for some n G Nq (lightning of odd order), then the solution 


ofLI'' + ^I = Ai{5^‘^'^+^^){t), 1(0) = /'(O) = 0, in £{R) i 


IS 


A 


I(t) = ^^-^(-l)V''i(5(2»-2fc))(o) cosa;t+ 


fc =0 


+ (- sinwt — i{H{t) sin cut))-!- 

n 

+ (_l)»^^2(n+l)^(^0)(0)cOSWt + ^(-l)V^(<5(2’^-2^))(t) 


k=0 


(5) IfV{t)=Ai{6(^^+^^){t)for some n G No (lightning of even order), then the solution 
of L I" + ^ I = ALi6^^^+^^){t), 1(0) =/'(0) =0, ini® 

Air 


I{t) = - [- ^(-l)V^-^t(5(2"+2-2^))(0) sina;t+ 


fc =0 


+ (—(-cos wt — l{H( t) cos oot)^ — 

n 

- (-l)"u;2’^+3t(Go)(0) sinwt + 

fc =0 

The derivation of the formula is similar to the one in the previous example and leave it 
to the reader. 


Our last example is outside the the scope of Theorem [3 and Theorem [51 

Example 8 (Lightning Rod under Lightning). Let A = | lns| = q{X,p), where A,^ = — InR^ ^ 
and A = q{A^) G M+. Notice that A is a positive infinitely large number, i.e. n < A for all n G N. 
ITe associate the initial value problem -^I” + + A^I = Ai{5'){t), 1(0) = /'(O) = 0, with the 
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LRC-electrical circle (lightning rod) with L = C = 1/A^ and R = 2/X and electromotive force 
V{t) = Ai{6){t) (under lightning of zero order with amplitude A). We prefer the form 
I" + 2X1' + A^/ = AX^ i{8'){t), /(O) = /'(O) = 0. For any test function (p € P(M) the initial 
value problem I" + 2X^p /' + A^ / = A^p A^ /(O) = I'{A) = 0, has a unique classical solution 

in £ (M) given by the formula 


Iip{t) — A^p Xp 


f g Kd gosa;(^(t — x) (/?(x) dx — — f e sinupft — x) p{x) dx— 

Jo Jo 


¥^(0) f . 

- - e ' sin ujpt 


00 , 




where oOp = A^^l — l/A^. This formula is obtained by Laplace transform method, but the 
reader is invited to verify its validity by direct differentiations and evaluations in the framework 
of £{M.). Notice that e~^A a moderate p-net since e~^‘e^ = {R^pY < {Rp)'"' for m = min{n G 
Z:n<i<n + 1} (Todorov h Vernaeve [13], Consequently {Ip{t)) is a moderate net. 
We also observe the the above formula implies that {Ip(t)) does not depend on the choice of the 
representative of the generalized number X modulo the ideal M{£{M.)^) (Todorov & Vernaeve [T3] . 
^f). Thus the generalized function I{t) = q{Ip{t)) € £^(]R) is the solution we are looking for. We 
present the final result in a “distributional-like” form: 


lit) = A 


- -e ' sinujt + Hit)e cos ujt — e - ) cos cot— 

Ul ^ iO ' 


— e 


—Xt 


XC A 

(S -) sinwi- H{t)e~^^ sinwt 

w ^ OJ 


where the generalized constants: ui,S,C in M \ M and the generalized functions: e sinwt, 
cosut , H(t)e~^^ sin ujt and Hit)e~^^ cos ut in^(R)\P'( 


are defined by 
0 


/ u pU 

iujpx) ipix) dx ), C = qi / cos ioj^x) </?(x) dx ), 

-CO j —CO 

e~^^ = s^ = qie~^'^^), sin ojt = qisinoj^t), cos uit = gicosLo^t), 


Hit)e sinujt = q 


(/ 


g-A^(t-a:) 


Hit)e cosiot = q 


— CO 

t 


Q-Kh-x) 


sinujpit — x) (fix) dx) and 
cos uipit — x) (fix) dx ), 


respectively. 
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